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^ ■ Abstract 

\ Erraticity analysis of multiparticle production data is introduced as a way of 

C*~) ' extracting the maximum amount of information on self-similar fluctuations. It 

is presented as the next logical step to take beyond the intermittency analysis. 
', An erraticity spectrum e(a) can be determined analogous to the multifractal 

<3\ ; spectrum f(a). An analytical example is presented to elucidate the method of 

r~| ■ analysis and the type of results that can be obtained. 



1 A historical overview 



Andrzej Bialas has played an influential role in the physics of hadronic and nuclear 
^ ■ collisions throughout his career. The work that he did with Robi Peschanski on in- 

termittency has dominated the attention of physicists working on multiparticle pro- 
duction in the last ten years. It is fitting at this point to review the significance of 
intermittency and ask where we can go from here. 

When many particles are produced in high-energy collisions, the very natural 
quantities to study theoretically and experimentally are averages, such as the mean 
multiplicity (n), the first few moments of the multiplicity distribution P n , and the 
rapidity distribution dn/dy. Indeed, those were the quantities investigated intensively 
in the beginning of the era of multiparticle production. 

Then as the collision energy was increased, the total rapidity range Y became 
large enough to permit meaningful partitioning of Y into smaller bins of various sizes 
5. It was found that the distributions P n (5) can be well fitted by negative binomial 
distributions with the normalized width increasing, as 8 decreases [0, |2|] . Thus began 
the interest in the study of multiplicity fluctuation as a function of the bin size. Such 
studies did not catch fire until the significance of intermittency, proposed by Bialas 
and Peschanski [H], was fully appreciated. 
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In particle physics intermittency refers to the power-law behavior of the normal- 
ized factorial moments F q , as the bin size is decreased. The observation of that 
behavior 0] therefore suggests that the mechanism for particle production has a self- 
similar property. It means that the occurrence of a large burst of particles in a small 
bin is rare, but possible if one waits long enough for such an event to take place. Any 
model that does not possess such intermittent features is thus ruled out. 

In recent years it was found that much of the intermittency phenomenon can 
be attributed to Bose-Einstein correlation among like-sign charged particles. The 
bunching of particles in small bins cannot be distinguished from the interference 
effect due to the coherent emission of same type particles from an extended source. 
While this is an important experimental finding, one should not let the BE correlation 
effect completely obscure the intermittency behavior, which is still seen in the unlike- 
sign charged-particle F 2 [||, |J . If it exists in the unlike-sign sector, then it must also 
exist in the like-sign sector, though small in comparison to the BE effect. An effect 
that is small is not necessarily unimportant. In this case it is our only clue to an 
important aspect of the dynamics of soft interaction: self-similarity. Hence, in my 
view the study of intermittency should go on. 

If the dynamics is self-similar, it is natural to ask about the multifractal properties 
of multiplicity distributions. The G moments were constructed to exhibit those prop- 
erties through the multifractal spectrum f{a) 0. The advantage is that the order 
q of the moments G q is continuous, and can be negative. That makes possible the 
determination of a which is a derivative in q, and facilitates the study of dips in ad- 
dition to spikes in the rapidity distribution. The disadvantage is that, unlike F q , the 
G q moments do not filter out the statistical fluctuations automatically and therefore 
require explicit elimination by "subtraction" H |J. When that is done, the dynam- 
ical Renyi dimensions D dyn can be determined. For the Monte Carlo code ECCO 
[ [Tcfl that simulates hadronic collisions with intermittency, it is found that D dyn jd is 
independent of the dimension d in which the G- moments analysis is done ||. A way 
to continue F q to noninteger values of q, while maintaining the attribute F q = 1 for 
all q when P n is Poissonian, has been devised Oj. Its application to real data has 
recently been attempted fl2 |. 

At this point there is a slow-down in the acquisition and analysis of multiparticle 
data of hadronic collisions. New methods of analysis have been proposed, notably 
by means of correlation integrals [H| and wavelets [fyj. They are more efficient and 
powerful than studying F q in discrete bins, and can extract more information on self- 
similiarity. The application of the wavelet analysis to real data has not yet been done, 
and the reward for such an improved analysis remains to be realized. 

While the progress in phenomenology is slow, one nevertheless can ask the the- 
oretical question: what is next? Is intermittency analysis the most that one can do 
to extract information of self-similarity of the particle production process? In the 
following section a suggestion is made to carry the study to yet another level where 
more information on fluctuations can be obtained. 
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2 Erraticity 



Let us examine in detail the normalized factorial moments F q . The horizontally 
averaged vertical moments are 

F M_ 1 f (n(n-l)---(n- g + l)) 3 . 

where (• • •) • is the (vertical) average over all events of the quantity bracketed at the 
jth bin, n being the multiplicity in that bin, and M is the total number of bins (e.g., 
M = Y/5 in the 1-dimensional case). If the space in which the partition into M bins 
is done is made to have a flat single-particle distribution by use of the cumulative 



variable |L5], jig] , one can also meaningfully define the vertically averaged horizontal 
moments 

m _ 1 £ (n(n-l)---(n-g + l)) e 

where (• • •) is now the (horizontal) average over all bins for the eth event, and M 
is the total number of events. It is clear that the two definitions are complementary 
and in most instances they behave the same way. In either case intermittency refers 
to the scaling behavior 

F q oc , (3) 

when M is increased in a fixed portion of the phase space, i.e., when the bin size 5 is 
decreased. 

The numerator of (|l|) and (|2|) are nonzero only when the bin multiplicity n is > q. 
Thus they pick out events and bins with large fluctuations, n ^> ( n )j e , when 5 is 
small, since (fi)- e oc 5. It is possible that, when q is large and 5 is small, one may 
have to wait for many non-contributing events to go by before finding a spike that con- 
tributes. That is why Bialas and Peschanski have coined the word intermittency for 
the phenomenon. The emphasis on bin multiplicity fluctuations marked a significant 
advance that intermittency generated in the subject of multiparticle production. 

However, intermittency as studied so far has not fully exhausted the character- 
ization of fluctuations that the system can exhibit. Let us focus on in (J^j to 
be definite. The summand is a quantity that characterizes the "spatial" fluctuations 
(in phase space or any other space) in an event. Since it plays a central role in the 
discussion to follow, let us denote it by F e so that 

q= We ' W 



' e 



We then see from (|]) that F^ 1 ' is an average of F q over all events. We know that F q 
fluctuates greatly from event to event. Those fluctuations are ignored by the study of 
F^ h \ so intermittency in F^ does not fully account for all the fluctuations that the 
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system exhibits. To capture the nature of those fluctuations and to find the associated 
scaling behavior constitute what can be called the erraticity analysis, which I now 
describe. 

It should be remarked that the problem to be addressed is not removed by up- 
grading to correlation integrals or wavelets. We shall use F^ as defined in ([J) as 
one possible, but simple, characterization of the spatial pattern of an event. Other 
descriptions can be chosen, and can be denoted by F^, used as a generic symbol. 
Indeed, F^ need not refer to multiparticle production. Any system that involves 
repeated samplings whose outcome can fluctuate from event to event can be inves- 
tigated in the erraticity analysis. To emphasize the generality of the method, let us 
simplify the symbol F^ to F e , when the order q is immaterial to the discussion. 

With F e describing the spatial pattern of an event, there should exist a distribution 
P(F) of F after many events. Let P(F) be normalized 

P(F)dF = 1 . (5) 
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Clearly, F^ in (§) is the average 

1 A/" „oo 

( F ) = JfT { F e = I FP(F)dF , (6) 

which conveys only a small piece of the information about P{F). Experimentally, the 
whole distribution P(F) should be determined. However, that may provide too much 
information, if the q and M dependences are fully explored. Thus a few moments of 
P(F) may be sufficient. Define the standard normalized moments 

where the averages are calculated as in (^). The order p here need not be an integer; 
in fact, it can even be less than 1, but may or may not be less than 0, depending on 
whether there are events that have F = 0. For p > 1, C p reflects the large F behavior 
of P(F), which is sensitive to the spikes in phase space. For p < 1, C p probes the 
low F behavior of P(F), which is influenced mainly by bins with low multiplicities, 
including empty bins. Thus knowing C p for < p < 2, say, reveals a great deal about 
the properties of P(F), all of which are not probed by intermittency. 

Collecting all the complicated properties of a complex system contributes only to 
a messy assemblage of facts. It is only when there is some simple, universal feature to 
be found in the assemblage that the phenomenological analysis becomes worthwhile. 
If there is self-similarity in the dynamics of particle production, we should search for 
power-law dependence on M. Eq. (3) exhibits one such behavior. Having generalized 
(F) to C p , it is natural for us to suggest the search for the scaling behavior of C p 

C p oc M^ v) oc . (8) 

While the behavior in @ has been referred to as intermittency, we shall refer to the 
behavior in (§) as erraticity. Since C p are the moments of P(F), they describe the 
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deviation of F e from the mean (F). Consequently, C p is sensitive to the erratic fluc- 
tuations of F e from event to event. Those fluctuations depend on the bin size because 
F e itself is a description of the spatial pattern that varies according to resolution. 
Thus if those fluctuations scale with bin size, then the erraticity exponent ip(jp) is an 
economical way of characterizing an aspect of the self-similar dynamics that has some 
order in its erratic fluctuations. 

Of particular interest is an index \x defined by 



a* = ±m 



(9) 

P =i 



It was shown in ||T7| , |18 that /i is related to the entropy in event space, and has been 
used to study chaotic behavior in branching processes in QCD. 

A dynamical system that has erratic fluctuations may or may not exhibit chaotic 
behavior in the technical sense of chaoticity in nonlinear dynamics. The generalization 
of the notion of chaos in classical trajectories to quantum systems where the degrees 
of freedom can increase with time is still under investigation. Whatever the outcome, 
the notion of erraticity is independent of it, and the results of erraticity analysis 
describe some features that are important in their own right. 

For a multifractal system one usually determines the multifractal spectrum f(a). 



It is related to a scaling exponent r(q) by a Legendre transform [OJ. In multiparticle 
production r{q) appears in the scaling law of the G moments @ 

G q {8) oc 5 T(q) (10) 

where q is a continuous variable. The exponent a is defined by 

<Mg) 



and the transform is 



f(a)=qa-r(q) . (12) 



It can be shown that for a multifractal set f(u) is always < a, and that the information 
dimension D\ corresponds to D\ — a\ — f(a±) at q — 1. 

Since C p is not the same as G q (their scaling laws (§) and (fL0|) having oppo- 
site behaviors in 5), our measure of erraticity does not have multifractal properties. 
Nevertheless, we still can define a spectrum e(a) by Legendre transform 

e(a) =pa - i/j(p) , (13) 



a„ = ^ . (14) 
dp 



The function e(a) exhibits certain properties of erraticity more directly than ip{p)- 
For example, we have oti = //, which is the only point where e(a) = a. For all other 
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values of a, one has e(a) > a. Since, by definition, C p = 1 at both p = and 1, we 
have i^(0) = = 0- Thus for any P(F) that becomes wider at smaller 5, as is 
always the case in particle production, it follows that ip(p) > for p > 1, ipijp) < 
for < p < 1, and 4>{p) > for p < 0, if C p exists. The resultant behavior of e(ct) is 
therefore that e(a) = at p = where cto < 0, and e(ct) > everywhere else, where 
a is calculable. An example of this behavior will be given in the following section. 

The values of and cti bear specific relationships to certain averages over P(F). 
Since (0) implies 

±C^j; iFl nUL) (JLJP(F) , (!5) 



it then follows from (||) and (14) that in the scaling region 



ao= l^M i{£nF) ~ £n{F)) ' (16) 

If the scaling laws are (F) oc M v and (£nF) oc ipinM, then we have 

"0 = ^-^ > (18) 

which is negative, except for unusual P{F). On the other hand, ( p!7[) suggests that if 
we define V = F/ (F), then 

(V£nV) = a 1 £nM (19) 

in the scaling region. The connection of ot\ with entropy should therefore not be 
surprising. 

In the foregoing we have suppressed the symbol q if F is the normalized factorial 
moment defined in (|j). For such moments to describe the spatial pattern, all the 
relevant quantities in (||) to (|19|) should be labeled with an index q, viz. , F q , C PtQ , 

i>q{p), fJ-q, e q( a ), a p,q, ^ <f V 



3 An analytic example 

To help make the discussion in the previous section more concrete and transparent, let 
us consider an example with analytic expressions. In real experiments or in computer 
simulations the event-to-event fluctuations of F may be so erratic that no simple 
formula can approximate P(F), let alone its dependence on q and M. However, 
there are general trends characteristic of multiparticle production that can be built 
in. Furthermore, Monte Carlo simulations of jet fragmentation in pQCD give definite 
shapes of P(F) that can serve as a very useful guide for the choice of analytic formulas. 



We shall rely on the results of |T8[ to generate specific expressions. 
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It should first be remarked that if the bin multiplicities in an event vary according 
to the Poisson distribution, then F q = 1 for all q. Fluctuations from F q = 1 have 
dynamical content, and usually (F q ) > 1, for q > 2. However, from the simulations 
in |Tj| we have seen that P{F q ) has its maximum at F q — 1 for all M. That is a 
condition that we shall impose. 

We adopt the gamma distribution for P(F): 

P(F) = AF a e~ bF , (20) 
and require that its peak be located at F = 1. With the normalization (|5|), it becomes 

W = f^T/"^ F ■ < 21 > 

where a is the only parameter, dependent on q and M. All the F distributions 
determined in |Tj| have the shapes of (|2l"D with large values of a. We adopt the 
following parametrization to introduce the q and M dependences: 

a = 500/(g InM) 2 , (22) 



which reproduces the general trend of the simulated results in ||18|| . 

In Fig. 1 are shown, as examples, the distributions P(F) for q = 2 and M = 
5, 50, 500. Clearly, for small bins (large M) there are large fluctuations of F from 
event to event. It is those fluctuations that are measured by the moments C p ^ q (M). 
For p = 2 the scaling behavior (§), as shown in Fig. 2, is satisfied at large M. Similar 
scaling behaviors are found for other values of p. Thus the exponents ip q (p) can be 
determined by the straightline fits in the scaling region. The results are shown by 
the dots in Fig. 3 for q = 2 and 3 and for some discrete values of p. Those dots are 
well fitted by the formula ip q (p) = J2t=i hp\ which is then used to determine a ps . It 
is clear from the general behavior of ip q (p), which is negative in the region < p < 1, 
that a , g < and a\ A > 0. Their numerical values are 

a 0i2 = -0.026, a ,3 = -0.047, 

a 1>2 = 0.024, ai, 3 = 0.04. 1 ' 

With analytical formulas for ip q {p) the spectrum e q (a) can be determined by use 
of (p~3|) . For the exponents ip q (p) shown in Fig. 3, the corresponding e q (a) are shown 
in Fig. 4. The straightline is for e q (a) = a. Thus where e q (a) curves touch the 
straightline are the values of and where e q (a) = give a ,g- The entropy indices 
are \i q = They get larger at higher values of q, which is a consequence of the 
fact that F q fluctuates more from event to event at higher q. 

Real data are not likely to be describable by simple formulas like (^) and (f2~2"l). 
However, erraticity analysis can be applied to the data directly, and curves for ip q (p) 
[using (0) and @] and e q (a) [using (pT3|) and (|14])1 can be determined, if scaling be- 



havior exists. That represents the "maximum" amount of information extractable 
from the horizontal and vertical fluctuations of the data that exhibit properties of 
self-similarity. Generally speaking, positive a describes the spikes of the spatial dis- 
tribution, while negative a describes the dips. 
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4 Conclusion 



Primitive averages are performed over both spatial fluctuations and event fluctua- 
tions, which have been referred to as horizontal and vertical averages, respectively. 
Intermittency probes the scaling properties of one of those fluctuations, and averaging 
over the other. Only one kind of moments are considered, viz., F q . Erraticity probes 
both types of fluctuations, and therefore double moments are needed: C p>q . 

Vertical fluctuations may be due to trivial reasons, such as impact parameter vari- 
ation from event to event. In heavy- ion collisions such variations should be controlled 
by Et cuts. For hadronic collisions, cuts in event multiplicity may restrict event 
fluctuations too much and unduly suppress the erraticity to be uncovered. Those 
fluctuations have dynamical as well as geometrical (i.e., impact-parameter related) 
origins and should be investigated fully. No geometrical fluctuation is present in 
e + e _ annihilations, so erraticity analysis explores the quantum fluctuations of parton 
branching for every fixed initial state specified by the energy. 

When the dynamics of particle production is known, erraticity analysis can then 
describe some aspects of that dynamics, such as the (possibly) chaotic behavior of 
perturbative QCD [|17], [18[. But the purpose of studying intermittency, and now 
erraticity, is to get phenomenological information from the data that can help us to 
learn more about the dynamics of particle production where the theory is inoperable. 
Specifically, it is for learning about the soft interaction. Plots of ip q (p) or e q (a) form 
the arena where theory and experiment should meet. Models of soft interaction that 
can reproduce the primitive averages may well reveal deficiencies when confronted 
with erraticity data. 

Heavy-ion data have so far not led to interesting results in intermittency study. 
Perhaps too much averaging has been done. Erraticity analysis may reveal more 
structure. 

Bose-Einstein correlation has temporarily detracted the study of intermittency. 
Focusing on unlike-sign charged particles in hadronic collisions and going deeper into 
erraticity may reveal features about the dynamics of soft interaction that may finally 
lead to the construction of a reliable model capable of meeting all experimental tests. 
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Figure Captions 

Fig. 1 Examples of the distribution P(F) given by Eqs. (21) and (22) for q = 2. 

Fig. 2 Scaling behaviors of C Piq (M) for p = 2 and q — 2, 3, 4. 

Fig. 3 Dependences of the erraticity exponents ip q (p) on p for q = 2, 3. 

Fig. 4 The erraticity spectrum e q (a) for g = 2, 3. 
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